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ABSTRACT 

Composite spectra of 85 proximate absorbers (log N(H i) > 20 crrT 2 and velocity differ- 
ence between the absorption and emission redshift, AV < 10,000 km s ) in the Sloan Digital 
Sky Survey are used to investigate the trends of metal line strengths with velocity separation 
from the QSO. We construct composites in 3 velocity bins: AV < 3000 km s , 3000 < AV 
< 6000 km s and AV > 6000 km s , with further sub-samples to investigate the metal line 
dependence on N(Hi) and QSO luminosity. Low (e.g. Sill and Fell) and high ionization (e.g. 
SilV and CIV) species alike have equivalent widths (EWs) that are larger by factors of 1 .5 - 3 
in the AV < 3000 km s composite, compared to the AV > 6000 km s _1 spectrum. The EWs 
show an even stronger dependence on AV if only the highest neutral hydrogen column density 
(log N(Hi) > 20.7) absorbers are considered. We conclude that PDLAs generally have higher 
metallicities than intervening absorbers, with the enhancement being a function of both AV 
and N(Hi). It is also found that absorbers near QSOs with lower rest-frame UV luminosities 
have significantly stronger metal lines. We speculate that absorbers near to high luminosity 
QSOs may have had their star formation prematurely quenched. Finally, we search for the sig- 
nature of dust reddening by the PDLAs, based on an analysis of the QSO continuum slopes 
relative to a control sample and determine a limit of E(B - V) < 0.014 for an SMC extinction 
curve. This work provides an empirical motivation for distinguishing between proximate and 
intervening DLAs, and establishes a connection between the QSO environment and galaxy 
properties at high redshifts. 

Key words: 



1 INTRODUCTION 

Damped Lyman alpha (DLA) systems with small velocity separa- 
tions (typically AV <3000 to 5000 km s ) from the systemic red- 
shift of the background QSO are traditionally excluded by most sta- 
tistical surveys (e.g. Lanzetta et al. 1991, 1995; Wolfe et al. 1995; 
Storrie-Lombardi et al. 1996; Ellison et al. 2001; Peroux et al. 
2001; Jorgenson et al. 2006; Ellison et al. 2008; Prochaska & Wolfe 
2009; Noterdaeme et al. 2009). The logic for the AV criterion was 
originally imposed to avoid absorbers that were not representative 
of the intervening population (e.g. those clustered around the QSO), 
or directly associated with the QSO host or its outflows. A fur- 
ther concern for small AV absorbers is that proximity to the intense 
QSO radiation field may alter the internal ionization balance of the 
galactic interstellar medium (ISM). In turn, this could lead to the re- 
quirement for complex ionization corrections when converting ob- 
served column densities into elemental abundances. Despite some 
notable exceptions (Prochaska et al. 2002a,b; Dessauges-Zavadsky 
et al. 2004, 2006; Milutinovic et al. 2010), ionization corrections in 
intervening DLAS can generally be ignored, even when the column 
density is relatively low. 



Although these proximity concerns have driven DLA selec- 
tion criteria for two decades, there is relatively little observational 
evidence to support the imposition of a velocity cut. Apart from a 
higher incidence by a factor of ~ 2-4 (Ellison et al. 2002; Russell, 
Ellison & Benn 2006; Prochaska, Hennawi & Herbert-Fort 2008), 
the only empirical evidence that the proximate DLAs (PDLAs) dif- 
fer from the intervening absorbers has been the higher incidence of 
Lyce emission detected in the DLA trough (M0ller & Warren 1993; 
M0ller, Warren & Fynbo 1998; Hennawi et al. 2009). The QSO's 
radiation does not manifest itself in gross differences in highly ion- 
ized ISM species. For example, Fox et al. (2009) found no evi- 
dence for a more frequent occurence of NV. The CIV and OVI col- 
umn densities of the proximate absorbers are also indistinguishable 
from the intervening population (Fox et al. 2007a, 2007b), although 
larger samples are needed to confirm this. It has also been con- 
cluded by Rix et al. (2007) (based on the few PDLAs with high 
resolution echelle spectra) that their chemical abundance patterns 
and overall metal enrichment are consistent with the intervening 
population. 

We have recently re-assessed the chemical properties of PD- 
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Figure 1. Distribution of absorber (absorption redshift, H i column density 
and velocity offset from the QSO systemic redshift) and QSO (luminosity) 
properties measured from the SDSS spectra for our sample of 85 PDLAs. 



LAs using a sample of 16 absorbers with AV < 3000 km s -1 using 
high resolution echelle data. In Ellison et al. (2010) we reported 
the first evidence that some PDLAs do exhibit both distinct metal 
enrichment and evidence for ionization by a hard radiation source. 
Specifically, it was found that although the proximate absorbers can 
indeed display a wide range of metallicities, at high N(H i) they are 
typically a factor of 3 more metal-rich than the intervening popu- 
lation. Although many of the proximate absorbers exhibit NV (in- 
cluding absorption with large velocity offsets), the most intriguing 
evidence for a hard radiation source came from non-solar abun- 
dances of S and Ar. Expansion of the PDLA data set is clearly re- 
quired to confirm these initial findings and to understand the case- 
by-case dependence on QSO proximity. 

In this paper, we take a complementary approach by lever- 
aging the large statistical power of the Sloan Digital Sky Survey 
(SDSS). Although the majority of DLA abundance studies are per- 
formed at resolutions R ~ 40,000, with high enough S/N it is pos- 
sible to detect even weak, unsaturated lines in spectra whose res- 
olution is an order of magnitude lower (e.g. Pettini et al. 1994). 
Our strategy here is to construct high S/N composite spectra with 
various criteria (e.g. AV, A'(Hi), QSO luminosity), in order to in- 
vestigate which quantities most significantly affect the metallicities 
of the PDLAs. 

Unless otherwise stated, we assume C1 M = 0.3, Q.a = 0.7 and 
H = 70 km s" 1 Mpc" 1 . 



2 SAMPLE SELECTION 



2.1 PDLA sample 

The largest compilation of PDLAs is the SDSS data release (DR) 5 
sample of Prochaska et al. (2008b, hereafter PHHF08). The parent 
sample used by PHHF08 to search for PDLAs excluded broad ab- 
sorption line (BAL) QSOs and required a S/N>4 over at least one 



region of 20 consecutive pixels. The resulting PDLA sample con- 
tains 108 absorbers with AV < 3000 km s" 1 and log JV(Hi) > 20.3. 
In addition to its size, a second motivation for the adoption of the 
PHHF08 sample is that PHHF08 have re-computed the systemic 
redshifts of those QSOs with proximate absorbers, accounting for 
the well-known systematic blueshifts that are incurred when us- 
ing rest frame UV lines such as Lycf and CIV (e.g. Gaskell 1982). 
Throughout this paper, we use the values of z a b s (determined from 
the centroid of the strongest low ion) and z cm reported by PHHF08 
to compute the values of AV. The revised values of z cm take into 
account systematic shifts between the observed emission lines and 
those that are generally considered to be reliable estimators of the 
systemic value (such as Balmer and oxygen recombination lines). 
Uncertainties on these systematic shifts are typically a few hundred 
km s _1 . In order to explore the impact of QSO proximity on DLAs, 
we extend the nominal A^H i) and A V criteria to include proximate 
absorbers with log N(Hi) > 20 and AV < 10,000 km s -1 . These 
cuts are motivated by the practical limitations of SDSS's resolu- 
tion and the availability of improved redshifts in PHHF08. There 
are 326 absorbers that fulfill the A^H i) and AV criteria. The sam- 
ple is reduced by requiring that the spectroscopic S/N ratio in the 
i-band (as reported in the SDSS fits header) exceeds 10. We also ex- 
clude one absorber towards a BAL QSO that was not flagged in the 
original filter by PHHF08 (J014049. 18-083942.5). The final sam- 
ple consists of 85 proximiate absorbers with log N(H i) > 20, A V < 
10,000 km s , S/N > 10 and improved estimates of z era (TableQ] 
Figure lATV Figure [TJ shows the distribution of QSO and absorber 
properties. 1 1 absorbers in the final sample have log A^H i) < 20.3. 
Although DLAs are more strictly absorbers with log A'(Hi) > 20.3, 
and 'proximate' absorbers usually considered to be within ~ 3000 
km s -1 of the QSO, for convenience, we will use the term 'PDLA' 
to refer to the 85 members of the sample studied in this paper. 

2.2 Control sample 

A handful of the strongest metal transitions can be seen in individ- 
ual SDSS spectra. In order to compare the equivalent width (EW) 
distributions of the Sill A 1526 line (which is not only strong but of- 
ten located redwards of the Lyo- forest) in the PDLAs and interven- 
ing DLA population, a sample of 85 intervening DLAs is assem- 
bled as a control sample. The control sample is drawn from the DR5 
DLA catalogue of Prochaska & Wolfe (2009), supplemented with 
a non-statistical sample of lower column density absorbers (recall 
that the PDLA sample extends to log A'(Hi) = 20). For each PDLA, 
the intervening DLA which is the closest simultaneous match in 
Zabs, A'(Hi) and QSO fibre i magnitude is determined. The fibre 
magnitude criterion is included in the matching process so that a 
similar spectral S/N is achieved. Once a match is made, the inter- 
vening DLA is not replaced and can not be considered as a match 
for a different PDLA. We also require that the Sill A 1526 line is 
redwards of the QSO's Lyo- line, i.e. z cm < 1.255 x (1 + z abs ) - 1. 
A Kolmogorov-Smirnov test (KS) is performed to ensure that the 
control sample properties are well matched to the PDLAs. The KS 
probabilities are 92% (z ab8 ), 83% (iV(Hi)) and 98% (i magnitude) 
indicating that there is no difference in the distribution of these 
properties between the DLAs and PDLAs. 

A second control sample is constructed in order to assess the 
extinction properties of the PDLAs via the comparison on the spec- 
tral slopes of QSOs with and without PDLAs (Section|5]l. The con- 
trol sample is drawn from the list of DR5 QSOs with z > 2, no 
PDLA and no strong BAL features, as determined from visual in- 
spection. We do not exclude QSOs with intervening DLAs because 
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QSO i'mag Plate MJD Fibre z em z abs logiV(Hi) AV(kms"') 



J012747.80+140543.2 


18.56 


0425 


51898 


278 


2.4877 


2.4416 


20.30 


3991.57 


J014214.74+002324.3 


18.08 


0401 


51788 


500 


3.3734 


3.3481 


20.40 


1733.61 


J023903.43-003850.8 


18.71 


0408 


51821 


134 


3.0782 


3.0185 


20.35 


4371.54 


J025518.58+004847.6 


18.87 


0410 


51816 


466 


3.9936 


3.9145 


21.40 


4771.41 


J073718.16+323631.5 


19.10 


0541 


51959 


151 


3.0175 


2.8926 


20.30 


9487.00 


J075901.28+284703.4 


19.02 


0859 


52317 


453 


2.8550 


2.8223 


21.05 


2532.00 


J080050.28+ 192058.9 


18.89 


1922 


53315 


152 


3.9533 


3.9465 


20.25 


448.52 


J08 1240.68+320808.6 


17.43 


0861 


52318 


333 


2.7045 


2.6259 


21.30 


6391.40 


J08 1256.05+563746.9 


18.80 


1872 


53386 


474 


3.3170 


3.2251 


20.00 


6425.95 


J081256.05+563746.9 


18.80 


1872 


53386 


474 


3.3170 


3.3387 


20.35 


-1310.54 


J081518.56+291 153.9 


19.33 


0930 


52618 


593 


4.2552 


4.2588 


20.50 


-205.44 


J082107.61+310751.2 


17.04 


0931 


52619 


491 


2.6193 


2.5347 


20.15 


7060.41 


J082531.88+263619.2 


18.93 


1267 


52932 


051 


2.5640 


2.5409 


20.45 


1933.80 


J082612.54+451355.7 


19.16 


0548 


51986 


167 


3.8170 


3.7068 


20.50 


6979.90 


J082638.59+5 15233.2 


17.02 


0442 


51882 


528 


2.8438 


2.8333 


20.80 


820.62 


J0835 10.92+065052.8 


18.37 


1297 


52963 


517 


3.9781 


3.9556 


20.35 


1359.00 


J083914.14+485125.7 


18.64 


0550 


51959 


460 


2.9673 


2.9692 


20.60 


-136.08 


J090017.61+490001.9 


18.72 


0765 


52254 


306 


3.2080 


3.2054 


20.95 


185.42 


J090033.49+421546.8 


16.78 


0831 


52294 


201 


3.2954 


3.2456 


20.30 


3477.11 


J090940.67+330347.6 


18.51 


1272 


52989 


022 


3.7835 


3.6583 


20.60 


7922.54 


J091210.35+054742.0 


18.09 


1194 


52703 


257 


3.2407 


3.1231 


20.35 


8398.34 


J091223.02+562128.5 


18.75 


0451 


51908 


312 


2.9816 


2.8894 


20.55 


7027.37 


J091548.90+ 302542.6 


18.85 


1936 


53330 


614 


3.1928 


3.0678 


20.45 


9091.90 


J092914.49+282529.1 


17.54 


1940 


53383 


441 


3.4044 


3.2629 


21.10 


9743.62 


J093019.58+423803.9 


18.86 


0870 


52325 


042 


3.7363 


3.6098 


20.20 


8119.55 


J094453.89+372840.2 


18.81 


1276 


53035 


172 


3.3397 


3.3356 


20.45 


262.81 


J095256.41+332939.0 


18.74 


1945 


53387 


032 


3.3962 


3.3468 


20.70 


3369.33 


J095744.46+330820.7 


18.15 


1948 


53388 


137 


4.2088 


4.1794 


20.45 


2039.99 


J095817.81+494618.3 


18.64 


1006 


52708 


474 


2.3555 


2.2909 


20.65 


5831.18 


J095937.1 1 + 131215.4 


17.43 


1744 


53055 


368 


4.0717 


3.9131 


20.15 


9540.28 


J100409.35+120256.5 


18.87 


1744 


53055 


109 


2.8804 


2.8006 


20.60 


6232.89 


J101725.88+61 1627.5 


18.13 


0771 


52370 


155 


2.8069 


2.7681 


20.60 


3073.19 


J10261 1.35+341459.9 


18.45 


1958 


53385 


449 


3.3991 


3.4141 


20.60 


-1021.19 


J102619.09+613628.8 


18.43 


0772 


52375 


227 


3.8442 


3.7853 


20.35 


3657.30 


J103403.87+380248.4 


18.15 


1998 


53433 


343 


3.5605 


3.5161 


20.70 


2948.23 


J104837.40-002813.6 


19.28 


0276 


51909 


313 


4.0070 


3.8880 


20.65 


7214.72 


J105123.03+354534.3 


18.54 


2090 


53463 


105 


4.8990 


4.8208 


20.40 


4003.30 


Jl 10855.47+120953.3 


18.62 


1604 


53078 


383 


3.6737 


3.5453 


20.75 


8302.28 


J111151.60+133235.9 


17.15 


1752 


53379 


242 


2.4328 


2.3820 


20.30 


4401.43 


J111611.73+411821.5 


18.04 


1439 


53003 


595 


2.9829 


2.9426 


20.20 


3012.79 


Jl 13002.34+1 15438.3 


18.39 


1606 


53055 


102 


3.3939 


3.3171 


20.20 


5226.94 


Jl 13008.19+535419.8 


17.60 


1014 


52707 


211 


3.0495 


2.9870 


20.25 


4575.67 


Jl 13130.41+604420.7 


17.60 


0776 


52319 


250 


2.9069 


2.8755 


20.50 


2428.55 


Jl 13354.89+022420.9 


18.92 


0513 


51989 


428 


3.9865 


3.9147 


20.65 


4332.45 


Jl 15526.34+351052.6 


18.85 


2035 


53436 


Oil 


2.8374 


2.7582 


20.70 


6247.59 


J120359.07+341 114.2 


18.58 


2099 


53469 


635 


3.7484 


3.6740 


20.60 


4718.10 


J120359.07+341 114.2 


18.58 


2099 


53469 


635 


3.7484 


3.6860 


20.60 


3917.08 


J121324.58+423538.5 


18.85 


1450 


53120 


455 


3.7663 


3.7646 


20.60 


56.65 


J 1 22040.23+092326. 8 


18.05 


1230 


52672 


066 


3.1462 


3.1322 


20.75 


978.39 


J123840.93+343703.3 


18.44 


2020 


53431 


490 


2.5718 


2.4712 


20.80 


8551.19 


J123937.56+343701.8 


18.41 


2020 


53431 


501 


2.4632 


2.4812 


21.10 


-1727.50 


J124138.32+461717.0 


18.57 


1455 


53089 


257 


2.7697 


2.6672 


20.70 


8300.70 


J124640.37+111302.9 


18.02 


1694 


53472 


085 


3.1475 


3.0975 


20.45 


3594.52 


J125659.47+301439.0 


18.71 


2011 


53499 


388 


2.9457 


2.8805 


20.70 


4967.33 


J125759.22-011130.2 


18.64 


0293 


51994 


127 


4.1120 


4.0208 


20.30 


5358.04 


J125832.14+290903.1 


18.95 


2011 


53499 


171 


3.4830 


3.3512 


20.55 


8928.92 


J130152.57-030729.3 


18.47 


0339 


51692 


277 


3.0706 


3.0510 


20.50 


1336.92 


J130426.15+120245.5 


17.93 


1696 


53116 


278 


2.9775 


2.9135 


20.55 


4835.33 


J130426.15+120245.5 


17.93 


1696 


53116 


278 


2.9775 


2.9289 


20.35 


3726.19 


J132235.12+584124.6 


18.63 


0959 


52411 


237 


2.8507 


2.8183 


20.50 


2511.27 



Table 1. 85 PDLAs used in this study. Magnitudes are in the Petrosian system and are corrected for Galactic extinction. 



4 Ellison et al. 



QSO 


i mag 


Plate 


MJD 


Fibre 


Zem 


Zabs 


log JV(Hi) 


AV(kms _1 ) 


T1 11724 69+11 5254 5 


18.54 


2109 


53468 


052 


3.1829 


3.1742 


21.30 


617.43 


TI'^SIOS 18-09^01X9 


18 70 


0914 


52721 


215 


2 4149 


2.3618 


20 30 


4647 59 


1115117 1 1+512825 5 


18.29 


1043 


52465 


492 


2.9153 


2.8346 


20.80 


6192.41 


1115801 97+014916 f) 


18.62 


0856 


52339 


196 


2.8937 


2.8535 


20.40 


31 13.30 


11 15842 92+652216 6 


18.70 


0498 


51984 


398 


3.1731 


3.0641 


20.30 


7576.96 


J ljJ7Ju.il T t tOj7i.J . 1 


18 44 


1671 


51446 


305 


2.3641 


2 2572 


20 40 


9647 59 


1141006 ^9+SQ9^19 ^ 

Jltl 7UU. JiT J "Z. J 1 > 


17 69 


0789 


52342 


368 


2 1189 


2 2467 


20 95 


6597 21 


1141950 54+082948 2 

J Itl 7Jv.J fc tTvOi7"0.i 


18.66 


1811 


53533 


198 


3.0344 


3.0504 


20.40 


_] 157.78 


1142446 29+054717 1 


17.91 


1827 


53531 


269 


2.3944 


2.3072 


20.40 


7751.18 


T1 441 27 66+475048 8 


18.80 


1674 


53464 


608 


3.1977 


3.2130 


21.05 


-1091.46 


IMSTJCJ ^4-009^7 S 

J 1 tJjt7.JJ T UUi-JJ / .J 


18 49 


0309 


51994 


423 


2 5417 


2 4416 


20 40 


8189 79 


1145429 65+004121 2 


19.42 


0309 


51994 


466 


2.6567 


2.5657 


20.15 


7667.94 


T1 50726 12+440649 2 


17.80 


1677 


53148 


424 


3.1133 


3.0642 


20.75 


1595 06 

-7 .J 7J .\JVJ 


T1 5241 1 15+410517 4 

J ULtl J.JJT7JUJJ / .1 


18.86 


1678 


53433 


614 


3.9196 


3.8809 


20.70 


2199.96 


J154153.46+315329.4 


17.53 


1581 


53149 


220 


2.5530 


2.4434 


20.85 


9388.11 


J155814.51+405337.0 


18.85 


1054 


52516 


602 


2.6403 


2.5521 


20.30 


7356.66 


J160413.97+395121.9 


17.95 


1055 


52761 


459 


3.1542 


3.1625 


21.75 


-656.45 


J164716.62+313254.4 


18.53 


1341 


52786 


025 


2.5173 


2.4936 


20.25 


2028.25 


J170353.98+362439.6 


18.91 


0820 


52438 


631 


2.4764 


2.4186 


20.25 


5046.92 


J173935.27+575201.7 


19.06 


0358 


51818 


605 


3.2081 


3.1357 


20.35 


5242.13 


J210025.03-064146.0 


18.21 


0637 


52174 


370 


3.1295 


3.0918 


21.05 


2729.34 


J21 1443.94-005532.7 


18.68 


0986 


52443 


185 


3.4245 


3.4420 


20.50 


-1211.24 


J223843.56+001647.9 


19.02 


0377 


52145 


560 


3.4499 


3.3654 


20.40 


5750.84 


J23 1543.56+ 145606.4 


18.58 


0745 


52258 


355 


3.3769 


3.2735 


20.30 


7149.85 


J232115.48+142131.5 


18.44 


0745 


52258 


232 


2.5539 


2.5729 


20.60 


-1607.98 



Table 1. Continued. 85 PDLAs used in this study. 

the PDLA sample will have the same statistical occurence of in- 
tervening absorbers. In order to reduce the statistical uncertainty in 
the properties of our control sample, multiple control spectra are 
matched to each PDLA spectrum. The best match for each QSO 
with a PDLA in absolute Petrosian i magnitude and z sm is taken 
from the list of non-PDLA QSOs (without replacement) resulting 
in a sample of 85 control QSOs. The KS probabilities that the i 
magnitudes and z cm of the PDLA and control QSOs are drawn from 
the same distribution are calculated. The process is repeated, with 
each iteration taking the next best simultaneous match in i and z cm 
until the KS probability of one of the two quantities drops below 
30%. This process yields 20 control QSOs for every PDLA QSO. 



3 COMPOSITE SPECTRA 



Before creating data stacks, the absorption redshift determined 
by PHHF08 was checked by fitting gaussians to three transitions 
in the SDSS spectra, where detected: OI ,11302, Sill ,11304 and 
Sill ,11526. The average of the metal line redshifts was adopted for 
spectral stacking. The mean offset between the values reported by 
PHHF08 and our mean metal line redshifts is 0.0002±0.001, where 
0.001 corresponds to 75 km s at z=3. The SDSS pixel scale is 69 
km s _1 per pixel, so the typical shifts are <K 1 pixel. As stated above, 
the A V values are still calculated based on the PHHF08 redshifts for 
easy reference and comparison to that work. 

All of the spectra are shifted to the rest-frame of the proxi- 
mate absorber as determined from the metal lines or using the orig- 
inal Lyo- redshift in 3 cases where no metals were detected. The 
spectra are normalized using continua that are fit using the Star- 
link software DIPSO. Composites are made by median combin- 
ing the individual normalized, rest-frame spectra. Using the median 



has three advantages over a mean combination. First, the median is 
more robust against spurious spectral features and a small number 
of absorbers with extreme absorption properties. Second, the un- 
absorbed continuum has a value very close to unity in the median 
composite. The mean composite is affected by other absorption (at 
different redshifts) resulting in a continuum level noticably below 
unity. Nonetheless, we have performed the analyses in this paper 
for both mean and median composites and the basic conclusions 
are unchanged. Figure [2] shows the stack of all 85 absorbers. The 
composite has a S/N ~ 150 at <2000 A, decreasing to -60 at 2500 
A. Note the different y-axis scales in the different panels to accentu- 
ate absorption feature detections. The S/N is high enough to detect 
relatively weak species such as Sill A 1808, Nill A 1741 and Znll A 
2026. 

In Figure [3] we divide the sample into three approximately 
equal-sized samples based on AV. The cuts are made at AV < 3000 
km s _I (29 absorbers), 3000 < AV < 6000 (27 absorbers) and 6000 
< AV < 10,000 km s -1 (29 absorbers). The median N(H i) of the 
absorbers in each of the velocity sub-samples is 20.6, 20.4 and 
20.5 respectively. PHHF08 quote typical 7V(H i) errors of 0. 15-0.20 
dex, indicating that the AV sub-sample differences are not signif- 
icant. We also compare each of the samples with a Kolmogorov- 
Smirnov test and confirm that there is no statistical distinction in the 
A'(H i) distributions. The velocity sub-samples are further divided 
by A'(H i). When dividing by A^H i) we do not aim to equalize the 
number of spectra contributing to each composite. Instead, we are 
motivated by the physical distinction that might arise as the ISM 
of the proximiate absorbers becomes more shielded. In Figures [5] 
and [6] we show the composite spectra of QSOs with log A'(Hi) < 
20.7 and log A'(Hi) ^ 20.7 absorbers respectively [for clarity, only 
the lowest (AV < 3000 km s~') and highest (AV > 6000 km s" 1 ) 
velocity sub-samples are shown]. 

Assuming that the N(H i) distributions of the three A V samples 
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Figure 2. Median-combined composite spectrum of 85 proximate absorbers with W(Hi) > 20 and AV < 10,000 km s '.The A:-axis is in the absorber rest-frame. 
Note the different y-axis scales in the different panels to accentuate absorption feature detections. 
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3.16 
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4 


41±7 


441+15 


102±20 


297±20 


315+25 


225±20 


1/48 


3.36 


AV > 6000 km s - \ 7V(Hi)>20.7 
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Table 2. Rest-frame equivalent widths for the full composite sample and three AV sub-samples. The penultimate column indicates the metallicity as a fraction 
of the solar value, according to the EW(SiII /il526)-metallicity relation of Prochaska et al. (2008a). The top line refers to the composite shown in Figure[2] 
The second section refers to the 3 AV composites shown in Figure[3] The third section refers to the low and high AV samples with low A'(Hi) shown in Figure 
[5] The fourth section refers to the low and high AV samples with high A'(Hi) shown in Figure|6] The fifth section refers to cuts in QSO luminosity shown in 
Figure[7] 



are comparable, for a given metallicity the EWs of the various metal 
species should also be comparable in the absence of any proxim- 
ity effects. From Figure [3] it is visually striking that the majority 
of absorption lines are strongest in the lowest AV composite. This 
is investigated quantitatively by measuring absorption line equiva- 
lent widths. EWs are determined using single gaussian fits within 
IRAF's splot task. For the CIV ,U1548, 1550 doublet the deblend- 
ing option was used since the lines have mild overlap in stacks with 
the highest CIV EWs. Deblending was also necessary for the OI A 
1302 transition in some of the composites (where blending with 
Sill A 1304 was evident). The rest-frame EW of several detected 
transitions are given in Tableland Figure [4] provides a further vi- 
sual representation of how the EWs vary as a function of AV. In 
most cases, errors in continuum placement dominate over statisti- 
cal error. The EW errors listed in Table[2]reflect the range of values 
obtained over multiple fit realizations rather than the formal error 
from the random noise. 



4 RESULTS FROM COMPOSITE SPECTRA 
4.1 Metallicity indicators 

Sill and Fell are commonly used as indicators of DLA metallicity 
due to the relatively large number of observable transitions with a 
range of oscillator strengths. The Fell A 1608 line is often unsat- 
urated in intervening DLAs. For comparable A'(Hi) distributions 
in the three sub-samples, and under the assumption that the line is 
unsaturated, the EW of Fell /11608 can therefore be used as a rela- 
tive metallicity indicator. From Table[2]we see that the EW of Fell 
/11608 is twice as large in the AV < 3000 km s -1 sample than the 
AV > 6000 km s -1 composite. For PDLAs with 3000 < AV < 6000 
km s"" 1 the strength of Fell /11608 is intermediate between the two 
samples, being ~ 20% larger than measured in the high AV com- 
posite. Figure [3] shows that the same trend is present for other Fell 
lines. The mean absorption redshifts of the composite spectra are 
given in Table [2] in order to indicate any effects that might be at- 
tributed to metallicity evolution. As a guide, Dessauges-Zavadsky 
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Figure 4. For three different absorption species, the rest-frame equivalent 
width is plotted as a function of AV. The top panel shows measurements 
from the composite spectra of the full PDLA sample (see also the second 
section in Table|2j. The middle and lower panels show the EWs measured 
in the composite of spectra of the low and high A'(Hi) sub-samples respec- 
tively (third and fourth sections fo Table[2}. 
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(green), AV > 6000 (red) km s . The panels show the strongest absorption features visible in the full composite shown in Figure|2] 



et al. (2009) find that the [Fe/H] in DLAs changes by 0.19 dex per 
unit redshift. 

Although almost always saturated, the EW of Sill /11526 can 
also be used as a metallicity diagnostic. The EWs of heavily satu- 
rated metal lines yield information about the velocity spread of the 
absorption (Ellison 2006; Ellison et al. 2008), which in turn is ex- 
pected to correlate with mass. Murphy et al. (2007) and Nestor et al. 
(2003) have both found that there is a correlation between Mgll EW 
and metallicity, and even unsaturated metal lines appear to exhibit a 
correlation between velocity spread and metallicity (e.g. Ledoux et 
al. 2006). Prochaska et al. (2008a) have calibrated the relationship 
between EW(SiII A 1526) and metallicity over two orders of magni- 
tude from 1/300 to 1/3 solar. For the values of EW (Sill A 1526) in 
Table|2]the calibration of Prochaska et al. (2008a) yields metallici- 
ties of 1/25, 1/37 and 1/43 of the solar value for the low, intermedi- 
ate and high AV samples respectively. If no ionization corrections 
are needed, the Fell and Sill lines indicate that the metallicity of 
PDLAs with AV < 3000 km s are typically twice as high as ab- 
sorbers with 6000 < AV < 10,000 km s -1 . 



We now consider the metal line strengths within the same ve- 
locity cuts, but now we additionally divide the sample by H i col- 
umn density (Figures[5]and[6]l. When only the log A^H i) < 20.7 ab- 
sorbers are included in the stack, the EW of Fell A 1608 is only 50% 
larger in the low A V sample than the high A V composite (compared 
to a factor of 2 for the full stack). The EW of Sill A 1526 indicates 
metallicities of 1/37 and 1/40 solar for the AV < 3000 km s" 1 and 
AV > 6000 km s _I respectively. The only species whose EW in- 
creases by more than 50% in the AV < 3000 km s -1 spectrum is 
CIV (see middle panel of Figure |4j. Conversely, there is a much 
larger difference in the strengths of all of the metal lines when only 
absorbers with log A'(Hi) > 20.7 are included in the composite. 
This result is visually striking in Figure [6] The larger difference 
between EWs in the low and high velocity stacks with log Af(Hi) 
>20.7 is caused mostly by a large increase (often a factor of two) in 
the EWs of transitions in the former. This is true of both low ioniza- 
tion species such as Sill A 1526 and high ionization species such as 
CIV A 1548. The EWs of both of these elements are known to scale 
with metallicity (Fox et al. 2007a; Prochaska et al. 2008a). Using 
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the calibration of Prochaska et al. (2008a), the Sill A 1526 EWs in- 
dicate metallicities of 1/12 and 1/45 solar for the low and high AV 
(N(Hi) > 20.7) composites respectively, even though the former 
has a slightly higher mean redshift (3.15 compared to 2.86, see Ta- 
ble[2]l. This is consistent with the finding of Ellison et al. (2010) that 
the PDLAs have a higher metallicity than the intervening popula- 
tion when the N(H i) is large. However, the difference between EWs 
in the low and high velocity stacks with log N(H i) >20.7 (Figure[6) 
is exagerated in the case of CIV and SilV, whose EWs in the high 
AF composites significantly decrease at higher A'(Hi). The lower 
column densities of CIV and SilV in the higher A'(H i) composite 
(for AV > 6000 km s~') is surprising. At these relative velocities we 
would consider these systems to be largely intervening and Fox et 
al. (2007a) found that there is no trend between N(CIV) and A'(Hi) 
in intervening DLAs. However, at CIV column densities above ~ 
14.5 the doublet becomes saturated and only a lower limit can be 
derived, a fairly common situation for DLAs (Fox et al. 2007a). 
Since Fox et al. (2007a) show that both the column density and ve- 
locity spread of CIV correlate with metallicity, the lower EW(CIV) 



at high A'(Hi) and high AV may be another manifestation of the 
paucity (at intervening redshifts) of high metallicity absorbers with 
high N(Hi) (e.g. Schaye 2001; Krumholz et al. 2009). 

4.2 Ionization indicators 

Ellison et al. (2010) suggested that sub-solar SH/Sill ratios may 
be an indicator of a hard ionizing source. SII is usually studied 
via the triplet at 1250, 1253 and 1259 A. In the SDSS data, the 
strongest line (A 1259) is blended with the Sill A 1260 line. The next 
strongest line (A 1253) is clearly detected in our 3 AV composites, 
but there may be some contribution from Lya absorption in the AV 
> 6000 km s _1 sample. We therefore only consider the SII A 1253 
absorption in the low and intermediate velocity samples. Table [2] 
shows that although the low AV composite has a slightly higher 
(by 9 mA) EW, the measurements are consistent within the errors. 
If the silicon abundance is ~ 50% higher in the low AV sample, 
the sulphur EWs indicate that the AV < 3000 km s _1 PDLAs have a 
SH/Sill ratio 50% lower than the 3000 < A V < 6000 km s -1 sample. 
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Figure 6. Median-combined composite spectrum of proximate absorbers with JV(Hi) > 20.7 split by AV: AV < 3000 km s 1 (blue), and AV > 6000 (red) 
km s . The panels show the strongest absorption features visible in the full composite shown in Figure[2] These are the same transitions shown in Figure[3] 



Similarly, the Sill EW in the A'(Hi) > 20.7 composite (third section 
of Table[2) has increased by a factor of 2.5 at low AV, whereas the 
SII EW has increased by only 20% (and is actually consistent with 
the 3000 < AF < 6000 km s~' sample within the errors). Higher 
SH/Sill ratios in PDLAs are consistent with the trend identified by 
Ellison et al. (2010). However, the relatively large uncertainties on 
the SII EWs mean that this result should be regarded as tentative. 

The species which exhibits the most striking difference be- 
tween the 3 AV composites is CIV. There is more than a factor of 
three difference between the EW of CIV A 1548 in the low and 
high AV composite spectra. This factor is even larger when only 
the log yV(H i) > 20.7 are included in the stack, where the EW in 
the low AV spectrum is almost a factor of 6 larger than the AV > 
6000 km s~ ] composite. The difference between the low and high 
AV samples is much smaller for the low Af(Hi) absorbers, only a 
factor of ~ 2.5. We see a similar behaviour for the SilV EWs: a fac- 
tor of ~ 3.5 difference between low and high AV composites when 
log N(H i) > 20.7, but only 50% difference for the low N(K i) stack. 
This is ostensibly surprising, as one might expect the lower N(H i) 



absorbers to be more susceptible to ionization effects, whereas high 
N(H i) absorbers might remain shielded even if they are relatively 
close to the QSO. 

NV is observed in the low A V composite but unfortunately can 
not be compared to the intermediate or high velocity spectra, be- 
cause of blending by the Lya forest. However, Ellison et al. (2010) 
have argued that strong NV in PDLAs may be at least partly due to 
higher metallicity, rather than (entirely) due to enhanced ionization 
(see also Fox et al. 2009). 

4.3 Dependence on QSO luminosity 

We investigate the effect of QSO luminosity on the composite spec- 
tra by dividing the sample of AV < 3000 km s -1 absorbers (of which 
there are 29) in half, based on the QSO flux at /1=1500 A in the 
QSO rest-frame. The median luminosity used to divide the QSOs 
into 'low' and 'high' luminosity samples is 4 x 10 43 ergs/s/A. The 
profiles of prominent absorption lines are shown for the low and 
high luminosity composites in Figure|7] The EWs of selected tran- 
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sitions are given in Table|2] The EWs of the low and high ionization 
species alike are larger in the lower QSO luminosity composite, de- 
spite an excellent agreement in the median N(H i) values of the two 
samples (20.6 and 20.5) and a KS test probability of 86% that they 
are drawn from the same distribution. The difference in the profiles 
in Figure |7]is therefore not due to different N(H i) distributions and 
is also unlikely to be due to an ionization effect for several reasons. 
First, OI is expected to be robust against ionization effects (e.g. 
Vladilo et al. 2001). Second, in the presence of either a hard (QSO) 
or a soft (stellar) ionizing spectrum, metallicities derived from Sill 
or Fell are expected to be over-estimated (Howk & Sembach 1999; 
Vladilo et al. 2001; Rix et al. 2007). Third, Ellison et al. (2010) 
show that [S/H] is under-estimated by at least a factor of three in 
more than half of the A'(Hi) < 20.7 PDLAs in their high resolu- 
tion study, yet the SII EWs in the low and high flux composites are 
consistent within the errors. 

York et al. (2006) and Vanden Berk et al. (2008) conducted a 
similar analysis to ours, constructing composite spectra from SDSS 
data of Mgll absorbers at intervening and proximate velocities, re- 
spectively. When their samples were divided by QSO magnitude 
(at i = 19.1 and M, = 26.5, for the two studies) the faint QSO 
composite was found to have stronger metal lines for both inter- 
vening and proximate absorbers. It was concluded that this effect 
was likely due to a selection bias since weaker Mgll systems are 
typically only detected towards brighter QSOs. The (mean) com- 
posite spectra of the bright and faint samples were therefore not 
constructed of absorbers with the same Mgll EW distribution. 

Such a bias is not likely to affect our selection of PDLAs. 
First, the identification of a PDLA relies on the broad, saturated 
Lya which is easily detected in spectra with much lower S/N than 
we require here (recall that PHHF08 required a S/N> 4, but our 
criterion is an i band S/N > 10). Second, we have checked that 
the iV(H i) distributions of the low and high luminosity samples are 
consistent. This can be seen visually in the top left panel of Fig- 
ureQand is confirmed with a KS test (described above). Therefore, 
although we might expect the overall S/N of our low luminosity 
composite to be lower, this does not affect the EW of the metal 
lines we measure in our median composite. Nonetheless, we check 
explicitly that the dependence of PDLA EWs on QSO luminosity is 
not due to a bias by constructing composites of AV >6000 km s -1 
absorbers again split by QSO luminosity. The EWs of the low and 
high ionization lines are now largely consistent with each other, as 
might be expected for these absorbers, which most studies would 
consider as intervening. However, a quantitative comparison of the 
low vs. high luminosities compsites of the AV> 6000 km s~' sam- 
ple is hampered by a non-negligible difference in the typical Hi 
column densities. Both the median and mean N(Hi) values differ 
by ~ 0.3 dex with a KS probability of only 6% that they are drawn 
from the same population. 

It is not clear why the PDLAs should have higher metal line 
EWs when they are within 3000 km s~' of lower luminosity QSOs, 
but it is consistent with an inverse relationship between QSO lumi- 
nosity and PDLA metallicity. One caveat is that although the N(H i) 
distributions are consistent between the low and high luminosity 
sample, the mean absorption redshift of the high luminosity sam- 
ple is Az = 0.5 higher than the low luminosity sample. Dessauges- 
Zavadsky et al. (2009) found that [Fe/H] changes by 0. 1 dex over a 
Az=0.5, whereas the logarithmic difference in metallicites inferred 
from the Sill EWs is 0.2 dex. Given the uncertainty of convert- 
ing Sill EW to metallicity, it is premature to conclude that there 
is a definite link between QSO luminosity and PDLA metallicity. 
However, the concept of interplay between an AGN and nearby 
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Figure 8. A zoom of the Lya absorption trough of the stacked spectrum 
(black solid line) of the 29 AV < 3000 km s" 1 PDLAs. The red dashed line 
shows a simulated stack where each of the 29 PDLAs has an emission line 
corresponding to Lt ya = 3 X 10 42 ergs/s superimposed at the centre of the 
Lya absorption trough. 

galactic ISM is not without precedent, either observationally or 
theoretically. Irwin et al. (1987) found that the radio jet axis of a 
nearby active galaxy was aligned with an HI tail in a close compan- 
ion. They suggested that the companion's HI gas was being ram- 
pressure stripped by the AGN jet. Alignment of the QSO outflow 
with the line-of-sight has been previously invoked to explain the 
higher incidence of transverse DLAs, relative to the proximate line 
of sight absorbers (Hennawi & Prochaska 2007). In a case study 
of one transverse absorber located only ~ 100 kpc from the QSO, 
Prochaska & Hennawi (2009) find extreme kinematics and consider 
the case for QSO-absorber interactions. The effect of HI stripping 
has been modelled by Fujita (2008) and confirmed as an effective 
mechanism for sufficiently high energy AGN. Such interactions 
could prematurely (or temporarily) shut-down the star formation 
leading to lower metallicities. Such a scenario might also explain 
why the lower N(H i) PDLAs have metallicities as low as the in- 
tervening population (e.g. Table[2] Ellison et al. 2010) whereas the 
high A'(H i) have been able to retain their gas and accrue their met- 
als. 

4.4 Lya emission in the DLA trough 

The detection of Lya emission superimposed on the DLA trough is 
a rare occurence, but one that is apparently more common in PD- 
LAs than DLAs (M0ller & Warren 1993; M0ller, Warren & Fynbo 
1998; Leibundgut & Robertson 1999; Hennawi et al. 2009). The 
luminosities associated with the Lya- emission are typically 5-20 
xlO 42 ergs/s and may show large velocity widths and spatial ex- 
tents of tens of kpc (Fynbo et al. 1999; Hennawi et al. 2009). 

One of the PDLAs in our sample has been recently identified 
by Hennawi et al. (2009) as one of the most luminous PDLA Lya 
emitters with a total Lya luminosity of ~ 4 x 10 43 ergs/s. Visual 
inspection of the rest of our sample (Figure lAll yields no other 
candidate Lya- emitter and there is no emission signal in the core of 
the composite spectrum. In order to determine an approximate limit 
for the emission luminosity in our combined spectrum we artifi- 
cally superimpose a gaussian emission line in the core of the PDLA 
trough in the observed SDSS spectrum for the 29 absorbers within 
3000 km s -1 of the QSO. The same Lya- luminosity is assumed for 
every PDLA, where the total integrated flux then depends only on 
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Figure 7. Median-combined composite spectrum of proximate absorbers with AV: < 3000 km s , split by QSO luminosity at A = 1500 A. The half with the 
highest luminosities (Z4500 > 4 X 10 43 ergs/s/A) are plotted in magenta and the half with the lowest luminosities (Z4500 < 4 X 10 43 ergs/s/A) are plotted in 
cyan. The panels show the strongest absorption features visible in the full composite shown in Figure[2] These are the same transitions shown in Figure[3] 



the luminosity distance to the absorber. The width of the emission 
line is set to cr= 100 km s _I in all cases. The spectra are then shifted 
to the rest-frame with a (1+z) 4 scaling to conserve the total flux and 
stacked to produce a new composite. Figure [8] shows the stacked 
spectrum of both the original 29 low AV PDLAs and an overlay 
of one of our simulated spectra with emission (in this case where 
the flux corresponds to a Lya luminosity of L Lya = 3 x 10 42 ers/s). 
Our tests demonstrate that the composite spectrum is sufficient to 
detect luminosities of a few x 10 42 ergs/s if such luminosities are 
commonplace in PDLAs. However, this limit is likely to be opti- 
mistic for 2 reasons. First, in previous detections of hya emission, 
the emission line often appears to be offset by a few hundred km s -1 
from the centre of the hya absorption trough. Stacking analyses are 
quite sensitive to such offsets (Ellison et al. 2000). Second, if the 
flux is spread out over many hundreds of km s _1 , as seen by Hen- 
nawi et al. (2009), the broader profile becomes more difficult to 
detect. In the absence of a well characterized distribution of line 



offsets and widths, it would be premature to attempt more sophisti- 
cated stacking simulations. 



5 DUST IN PDLAS 

Over the past two decades, there has been considerable debate in 
the literature regarding the amount of dust extinction associated 
with DLAs. Although these studies may disagree in the details, 
there is now a broad concensus in the contemporary literature that 
the typical dust reddening is small. Some examples of recent works 
on this topic and detections or limits on the reddening include Mur- 
phy & Liske (2004) [ E(B-V) < 0.01], Ellison, Hall & Lira (2005) 
[ E(B-V) < 0.04], Vladilo et al. (2008) [ E(B-V) = 0.006] and 
Frank & Peroux (2010) who find a best fit E(B-V) = -0.002, i.e. 
consistent with no reddening. 

In this section, we investigate whether the PDLAs have extinc- 
tion properties that are consistent with the intervening population, 



12 Ellison et al. 



since there are several a priori reasons to expect that they may dif- 
fer. As we have demonstrated above, the PDLAs (particularly those 
at low AV with high A'(Hi)) have higher metallicities, which appear 
to go hand-in-hand with higher gas phase depletion indicators (e.g. 
Prochaska & Wolfe 2002). Moreover, Kaplan et al. (2010) have 
shown that their sample of metal-strong DLAs (a sample which 
overlaps considerably with our PDLAs, as we discuss in the next 
section) have systematically redder spectra than a sample of control 
(non-metal strong) DLAs. Finally, proximate Mgll absorbers at z ~ 
1-2 have reddening values that are approximately twice that of the 
intervening Mgll systems: E(B-V) = 0.02 and 0.01 respectively 
(Vanden Berk et al. 2008; York et al. 2006). 

There are three main techniques for determining the redden- 
ing in a sample of DLAs. The first uses broad band colour infor- 
mation, usually in the optical (e.g. Khare et al. 2004; Vladilo et al. 
2006, 2008), but if IR photometry is also available, a more strin- 
gent limit can be obtained (e.g. Ellison et al. 2005). If spectral data 
are available, a more detailed extinction analysis is possible and 
there are two main techniques that have been applied. The first in- 
volves fitting a power law to unabsorbed regions of the QSO con- 
tinuum, avoiding the emission features (e.g. Pei, Fall & Bechtold 
1991; Murphy & Liske 2004; Kaplan et al. 2010). The second in- 
volves making composite spectra of an absorber and non-absorber 
sample, where the ratio reveals the average extinction curve of the 
former (e.g. York et al. 2006; Wild et al. 2006; Vanden Berk et 
al. 2008; Frank & Peroux 2010). We have experimented with both 
techniques. The composite method was found to be fundamentally 
limited (in our relatively small sample) by our ability to average 
out the range of underlying QSO properties. Although this tech- 
nique has been used successfully on smaller samples than ours 
(e.g. Wild et al. 2006), the redshift range of our absorbers is such 
that there are many emission lines in our rest-frame spectrum. As 
noted by Frank & Peroux (2010), this leads to residuals and struc- 
ture in the region 1300-1700 A. Despite various attempts at match- 
ing parameters, we were unable to overcome this limitation and 
were always dominated by the underlying variation in emission line 
properties. This may be mitigated in larger samples (e.g. Frank & 
Peroux 2010) where residuals lead to an apparent bluing of only 
E(B-V)~ -0.002. 

Instead, we adopted the approach of fitting emission- and 
absorption-free regions of the QSO continua with a power law and 
examining the difference in spectral indices of the PDLAs and a 
control sample of no-PDLA QSOs. The technique is the same as we 
have applied previously in Kaplan et al. (2010). In brief, the median 
flux of up to seven emission-free regions (1312-1328, 1345-1365, 
1430-1475, 1680-1700, 2020-2040, 2150-2170, and 2190-2250 
A) is calculated from the rest-frame QSO spectrum. The number 
of regions depends on the QSO redshift: Higher redshift QSOs will 
have few spectral regions to fit, as the redder regions are shifted 
beyond the SDSS spectral coverage. A power law of the form 
/(/t) oc A~ a is fit to these median values where reddening due to 
dust will lead to smaller values of a. An example of the fitting pro- 
cedure is shown in Figure [9] Three QSOs have 2 PDLAs in their 
line of sight. These were excluded from the analysis, leaving a sam- 
ple of 79 PDLAs. Fitting is also performed for the 20 control QSOs 
(Section|2j2} matched to each of the 79 PDLAs (i.e. a total of 1580 
control spectra). The use of control spectra matched in redshift and 
magnitude will help circumvent any bias from the dependence of 
continuum slope on luminosity or the number of wavelength win- 
dows fitted. We have additionally checked that there is no overall 
trend of a with redshift in the control sample. 

The top panel of Figure [10] shows the fractional distribution 
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Figure 9. Example of the spectral fitting procedure. A power law (/(/I) cc 
A~ a ) fit is made to the median value (small green crosses) in each of up to 
seven continuum regions (vertical blue shaded bands). The solid black line 
shows the fit and the yellow region shows the uncertainty as a function of 
rest-frame wavelength. 





KS probability 31.45 % 




Measured PDLA indices 


6 


I Median PDLA a=1.35 

Median control a=1.35 










CO 

6 












o 




1 






- 


d 






















. . . 1 


o 


1 








2 ; 



Index (a) 





KS probability 0.196 % 




Measured PDLA indices — 0.09 ' 


<* 

d 


1 Median PDLA a=1.26 

Median control a— 1.35 










CO 

o 








— - 




d 






i ~ 
i 




- 


d 










: 















-10 12 3 



Index (a) 

Figure 10. Fractional distribution of power law indices for the PDLA sam- 
ple (red solid histogram) and matched control QSO sample (black dashed 
histogram). There are 20 control QSOs for every QSO with a PDLA. The 
upper panel shows the distribution of measured PDLA indices. The lower 
panel shows the PDLA indices decreased by a value of 0.09 which is the 
offset required to yield a KS value less than 0.3 % (3 cr rejection of the null 
hypothesis). The control histrogram is the same in both panels. 
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of power law indices in the PDLA and control sample. A KS test 
indicates that there is no statistical difference between the distribu- 
tions and the two samples have the same median value. In order to 
determine a limit for the reddening in our sample, we estimate the 
minimum change in the indices of the PDLA sample that would 
have resulted in a KS test probability of 0.3% that it is drawn from 
the same distribution as the control (corresponding to a 3 cr cer- 
tainty). Small incremental increases are made to the observed in- 
dices of the QSOs with PDLAs and the KS probability (relative to 
the control) is re-computed until the value reaches 0.3%. Consid- 
ering all 79 PDLAs, decreasing all the indices by a single value of 
0.09 yields a KS probability below 0.3%. The new distribution of 
PDLA indices is shown in the lower panel of Figure [Tol 

The sample of 79 PDLAs includes absorbers with AV values 
as large as 10,000 km s -1 , At such large velocity separations the 
absorbers are likely to be dominated by the intervening population. 
However, we find that comparing only the lower AV subsamples 
with their control still yields no significant difference in the dis- 
tribution of a. There is also no difference in the continuum slopes 
when the sample is split by A'(Hi). We repeat the above procedure 
for determining the minimum difference in a that would result in a 
significant KS result. For the AV < 6000 km s _I absorbers (51 PD- 
LAs) we again find an index decrease of 0.09 before the KS proba- 
bility drops below 0.3%. For the AV < 3000 km s" 1 (28 absorbers) 
the test is less sensitive and an index change of 0.19 is required 
before the distribution is significantly different from the control. 
Using the smaller data set therefore provides the more stringent 
constraint on reddening. 

To convert the index change that results in a significant KS 
probability into a measure of dust extinction, we define the differ- 
ence between the control and PDLA indices as 6a = a c - ap. As 
described in Kaplan et al. (2010), we can write the colour excess at 
two wavelengths as a function of 6a: 



-2.5.log h 



(1) 



The values of Ai and A 2 are set to be 1300 and 2000 A respec- 
tively, as values that are typical of the coverage of our spectra. In 
order to convert this colour excess at A\ and A 2 into the more stan- 
dard notation of E(B - V), we adopt the SMC extinction curve of 
Pei (1992). We define %(X) = A A /A V , so thaQ 



E Xl -x 2 =A M -A, h _ =A v (£ t] -£t 2 ). 



(2) 



Combining these equations for a 6a = 0.09 and 0. 19 we derive 
an A v = 0.020 and 0.041 respectively. To convert to an E(B- V) we 
use the standard definition R v = A V /E(B- V) with R v = 2.93. For a 
6a = 0.09 and 0.19 we arrive at upper limits to the colour excesses 
of E(fi - V) = 0.007 (AV< 10,000 km s" 1 and AV< 6000 km s" 1 
samples) and 0.014 (AV< 3000 km s -1 sample) respectively. 

An alternative technique for identifying differences in the 
PDLA and control sample is to calculate the difference between the 
a of the PDLA spectrum and the mean of its 20 control galaxies. 
These 'a offsets' are plotted as a function of AV in Figure [PT1 and 
should be zero if the QSO with the PDLA has a continuum slope 
that is the same as the mean of its controls. Negative values of the 
a offset indicate that the PDLAs are redder than their controls, on 
average. Figure [TT] appears to show that QSOs with PDLAs whose 
AV < 2000 km s are bluer than their matched controls. At larger 

' Pei (1992) defines ^(A) = A^/Ag, so we have included a factor of 1.33 to 
the values of £ derived from the Pei extinction curve. 
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Figure 11. a offset is defined as the index of the PDLA QSO spectrum 
minus the mean value of the index of its 20 control QSOs. The mean a offset 
and lcr error bars are plotted in bins of AV. Based on a random sampling of 
the full a offset distribution, the high value seen in the AV < 2000 km s -1 
bin has a 5% chance of occurring by chance in a sample of this size. 



AV there is no significant (> 3tr) offset. Even though we take the 
median flux in the continuum windows to account for the effect of 
absorption features, one may be concerned that indices may be bi- 
ased in certain velocity intervals as strong absorption lines move 
through. In order to make the PDLA spectrum bluer, additional ab- 
sorption at the red end of the spectrum would be required. A review 
of Figure|2]shows that the only strong lines beyond 1700 A are the 
Fell lines at Aq > 2300 A and are therefore too red to ever enter 
any of our continuum windows. Intrinsic dependences of the index 
on QSO redshift and magnitude are accounted for in the matching 
procedure. Drawing samples at random from the full a offset dis- 
tribution indicates that the positive offset seen at AV < 2000 km s _1 
would be expected 5% of the time for bins of this size. Larger sam- 
ples are therefore required to determine whether this result is due 
to small number statistics. 



6 DISCUSSION 

A similar analysis to the one presented here has been performed 
for 415 proximate Mgll absorbers (EW > 0.3 A) within AV < 3000 
km s" 1 of the QSO by Vanden Berk et al. (2008). They find that the 
low ionization lines (Mgll and A1II) have similar equivalent widths 
in the proximate and intervening composites, whereas the high ion- 
ization lines (CIV and SilV) are stronger in the AV < 3000 km s -1 
sample. Vanden Berk et al (2008) conclude that this is caused by 
higher ionization in the proximate Mgll systems. Although we also 
find that the CIV and SilV EWs are 2-3 larger in the PDLAs at AV 

< 3000 km s -1 , compared to the AV > 6000 km s -1 sample, the low 
ionization lines such as Sill and Fell are also up to a factor of 2 
stronger. The relative strengths of low and high ionization lines is 
probably related to the ,/V(Hi) of the sample. The Mgll absorbers 
in Vanden Berk et al. (2008) are likely to have H i column densi- 
ties that extend down to log A'(Hi) ~ 17.5. Their composites will 
therefore be skewed towards lower A'(Hi) values than ours. When 
we consider the lower half of our A'(Hi) distribution (log A'(Hi) 

< 20.7) we also find that the EWs of the low ions become more 
consistent in the low and high AV samples (Figures[4]and[5]l. The 
clearest distinction in the metallicities of AV < 3000 km s relative 
to larger velocity separations is also seen when the N(H i) is large 
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(see also Ellison et al. 2010). This may again explain why Vanden 
Berk et al. (2008) did not find evidence for enhanced metallicities 
in their proximate Mgll sample, whereas our PDLAs show system- 
atically larger low ionization metal line EWs at small AV. 

Figure [3] and Table [2] show that the majority of the transi- 
tions in the 3000 < AV < 6000 km s" 1 composite exhibit EWs 
that are intermediate between the low and high velocity compos- 
ites. This indicates that the 3000 - 6000 km s -1 absorbers may not 
yet have reached the same typical values as intervening absorbers. 
We can attempt to achieve some finer AV resolution by utilizing 
the Sill A 1526 line which is detected in 72/85 of the individual 
SDSS PDLA spectra. In cases of a non-detection, the rest-frame 
3 <t limit was calculated from the rest-frame spectrum using the 
equation 3cr = 3x ™" M ■ The FWHM at 1526 A is taken to be 0.85 
A (based on the nominal SDSS resolution of R = A/AA = 1800). 
The EW of this (usually saturated) Sill line correlates with metal- 
licity (Prochaska et al. 2008a) so can be used as an approximate 
metallicity metric. In Figure [T2] we show the median Sill A 1526 
EW in bins of AV (middle panel) and the corresponding metallicity 
calibrated by the Prochaska et al. (2008a) relation (bottom panel). 
The enhanced metallicity is seen most convincingly in the lowest 
velocity bin (AV< km s~') and is due to a combination of a lack 
of low EWs, as well as values that range up to 1 A in equivalent 
width. The results earlier in this paper would further lead us to ex- 
pect that this will depend on the N(H i) of the individual absorbers, 
but we do not have the numbers to test both the AV and A'(Hi) 
distributions separately in this way. 

Although the median Sill EW is only enhanced at AV < 
km s -1 , Figure [12] shows a notable population of high EW Sill ab- 
sorbers at larger velocity separations, particularly around 3000 < 
AV < 5000 km s -1 . These are likely to be the absorbers that lead 
to the generally larger EWs measured in the intermediate velocity 
composite (3000 < AV < 6000 km s" 1 ). Although Sill EW is not a 
very accurate estimator of metallicity, the detection of Sill A 1808 
(with EWs> 150 mA) in approximately half of the 3000 < AV < 
5000 km s -1 PDLAs with Sill EW > 900 mA is strongly sugges- 
tive of high metallicities. It may be surprising that QSO proxim- 
ity is responsible for enhanced metallicities at AV > 3000 km s -1 
(which corresponds to a Hubble flow distance of ~ 10 proper Mpc 
at z = 3). Although the redshifts derived from the SDSS spectra by 
PHHF08 are improved from those determined from Lyor and CIV 
emission lines alone, they are likely still only statistically accurate 
to within a few hundred km s -1 (Shen et al. 2007). However, the 
presence of a population of PDLAs with -2000 < A V < indicates 
that a combination of non-Hubble flow velocities and redshift er- 
rors can account of ~ 2000 km s in velocity range. Individual IR 
spectra of these QSOs are required in order to determine accurate 
relative velocities. The uncertain nature of these intermediate AV 
absorbers means that, even with a sample of 85 SDSS spectra, it 
is not possible to determine what the correct cut-off should be for 
studies that aim to sample the truly intervening population. Our re- 
sults indicate that the largest effect occurs at velocities well below 
the cut-offs typically used in DLAs surveys (~ 3000 km s -1 ). How- 
ever, a more conservative approach (until a larger statistical study 
can be undertaken) would be to increase the AV cut-off from 3000 
to 6000 km s -1 . 

Although the majority of DLAs are metal-poor, a few rare ex- 
amples of high metallicities have been reported in the literature, al- 
though mostly in the lower N(H i) sub-DLAs (e.g. Prochaska et al. 
2006; Peroux et al. 2008; Meiring et al. 2008; Dessauges-Zavadsky, 
Ellison & Murphy 2009). The very large samples of QSOs assem- 
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Figure 12. Top panel: Distribution of Sill A 1526 EWs for the PDLAs (solid 
red line ) and matched control sample of intervening DLAs (dashed black 
line). Middle panel: Rest-frame EWs of the Sill A 1526 line as a function of 
velocity offset from the QSO. Limits are 3<x. The open circles show individ- 
ual PDLA measurements and the solid circles are median values in velocity 
bins of width 2000 km s -1 . The median is selected so that the small number 
of Sill upper limits can be considered as detections without affecting the 
binned value. The horizontal dashed line is the median value of Sill A 1526 
EWs of the intervening control sample. Lower panel: Sill A 1526 EWs have 
been converted to metallicities (relative to the Sun) using the calibration of 
Prochaska et al. (2008a): log Z/Z Q = 1.41 log(W) - 0.92, where W is the 
Sill EW in A. The horizontal dashed line is the median metallicity calcu- 
lated from the Sill A 1526 EWs of the intervening control sample. Errors 
are lo\ 



bled by the SDSS has also permitted the identification of unusu- 
ally high metallicity systems. Herbert-Fort et al. (2006) identified 
absorbers with large metal line equivalent widths from the SDSS 
DR3 spectra. These metal-strong (MS) DLAs are of particular in- 
terest for detecting rare elements, such as boron, germanium and 
cobalt (Ellison, Ryan & Prochaska 2001; Prochaska et al. 2003). It 
is estimated that this metal-strong population represents only ~ 5% 
of the DLA population. Kaplan et al. (2010) have determined abun- 
dances for a sample of 16 MS DLAs finding a median metallicity 
Z ~ 1 /5Z at z ~ 2, a factor of around 4 higher than their control 
sample of non-MS DLAs. We have compared our sample of 85 PD- 
LAs (TableQ} to the list of MS DLAs in Herbert-Fort et al. (2006). 
1 1 PDLAs (13%) are present in the compilation of Herbert-Fort et 
al. (2006), see Tabled However, only 45/85 of the QSOs in Table 
[TJwere present in the DR3 (recall that our sample is drawn from 
the DR5). Therefore, whereas only ~ 5% of intervening DLAs are 
metal-strong, the percentage rises to 1 1/45=24% for PDLAs. 

One of the PDLAs in our sample is the proto-typical MS DLA 
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QSO z em z abs logAT(Hi) AVfkms- 1 ) 



J075901. 28+284703.4 


2.8550 


2.8223 


21.05 


2532.00 


J081240.68+3208O8.6 


2.7045 


2.6259 


21.30 


6391.40 


J095817.81+494618.3 


2.3555 


2.2909 


20.65 


5831.18 


J101725.88+61 1627.5 


2.8069 


2.7681 


20.60 


3073.19 


J102619.09+613628.8 


3.8442 


3.7853 


20.35 


3657.30 


J113008.19+535419.8 


3.0495 


2.9870 


20.25 


4575.67 


J122040.23+092326.8 


3.1462 


3.1322 


20.75 


978.39 


J135305. 18-025018.2 


2.4149 


2.3618 


20.30 


4647.59 


J160413.97+395121.9 


3.1542 


3.1625 


21.75 


-656.45 


J210025.03-064146.0 


3.1295 


3.0918 


21.05 


2729.34 


J2321 15.48+142131.5 


2.5539 


2.5729 


20.60 


-1607.98 



Table 3. PDLAs identified as 'metal-strong' in the DR3 catalogue of 
Herbert-Fort et al. (2006). 



J08 1240.6+32080^1 In a case study of this absorber, Prochaska 
et al. (2003) used data obtained with HIRES to make the first de- 
tections of several elements never before detected outside the lo- 
cal group. J08 1240.6+320808 is notable because it shows Sill A 
1808 (EW ~ 260 mA) absorption even in the SDSS spectrum. As 
described above, Sill A 1808 is observed in several of the individ- 
ual SDSS spectra, including J075901+284703, Jl 13008+535419, 
J135305-025018, J160413.97+395121 and J210025.03-064146 
which appear in the MS DLA catalogue of Herbert-Fort et 
al. (2006). To demonstrate the large EWs of these metal- 
strong PDLAs, in Figure [T3] we show the SDSS spectrum of 
J160413.97+395121 compared with that of J081240.6+320808. 
J 1 604 1 3 . 97+395 121 has been subsequently studied at high reso- 
lution by Ellison et al. (2010), who confirm a relatively high metal- 
licity ~ 1/16 Z Q . J210025.03-064146.0 was observed with ESI by 
Herbert-Fort et al. (2006) and measured to have an even higher 
metallicity: Z ~ 1/5 Z Q . 



Although the PDLAs and MS DLAs may both preferentially 
sample the high metallicity end of the DLA distribution, they differ 
markedly in their extinction properties. Whereas the general DLA 
population causes very little reddening to the background QSO 
(Murphy & Liske 2004; Ellison et al. 2005; Vladilo et al. 2008; 
Frank & Peroux 2010), Kaplan et al (2010) find significantly flatter 
power law slopes in QSOs behind MS DLAs. Using the method- 
ology described in Section|5]and adopting our choice of SMC ex- 
tinction curve and R v , their report value of 6a = 0.29 corresponds 
to A v = 0.063 and E(6 - V) = 0.022. In contrast, our most strin- 
gent limit (derived for the AV < 3000 km s sample) gives A v < 
0.041 and E(B - V) < 0.014 for the PDLAs. The reddening by 
PDLAs is also less than for associated Mgll systems at z ~ 1-2 
[E(B-V) = 0.02, Vanden Berk et al. 2008], although this may be 
at least in part due to the redshift evolution of reddening in Mgll 
absorbers (Menard et al. 2008). However, our detection limit is not 
sufficiently sensitive to determine whether or not the PDLAs are 
more or less dusty than their intervening cousins. 

We have presented evidence in this paper and in Ellison et al. 
(2010) that some PDLAs (those with small AV and large A'(Hi)) 
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Figure 13. Selected regions of the SDSS spectrum of 
J160413.97+395121.9 (black) with the SDSS spectrum of the proto- 
typical metal-strong DLA J081240.68+320808.6 (grey) overlaid. The top 
panel includes detections of SII AA 1250, 1259, Sill A 1260, OI A 1302 and 
Sill A 1304. The middle panel includes detections of Sill il526, CIV AA 
1548, 1550 and Fell A 1608. The lower panel covers Sill A 1808 and Al III 
AA 1854, 1862. 

appear more metal-rich than the intervening population. However, 
the sample of PDLAs in Ellison et al. (2010) includes an absorber 
with one of the lowest metallicities measured to date (J0140-0839, 
[0/H]=-2.72). In Section@3]we have shown that for PDLAs with 
AV < 3000 km s _1 , the composite SDSS spectra show weaker metal 
lines when the QSO luminosity is higher. J0140-0839 (AV= 1250 
km s~') is included in the PDLA sample of PHHF08, but excluded 
by us because the QSO shows mild BAL features. Nonetheless 
it is intriguing that J0140-0839 has a relatively large luminosity: 
^4500 = 1.3 x 10 44 ergs/s/A, a value exceeded by only two of the 
QSOs in our sample. There are 4 QSOs in our sample with L 1500 > 
1 x 10 44 ergs/s/A: J082638.59+5 15233.2, J090033.49+421546.8, 
J092914.49+282529.1 and J095937. 11 + 131215.4). The latter two 
have very large relative velocities from the QSO (AV > 9500 
km s- 1 ). We inspected the metal lines of J082638.59+515233.2, 
J090033.49+421546.8 (AV = 820 and 3477 km s _I respectively) 
in the SDSS spectra and did not find them to be abnormally low. 
High luminosities and moderate relative velocities do therefore not 
necessarily lead to low metallicities. 



2 We emphasize again that for convenience, in this paper we are using the 
term 'PDLA' for absorbers with log N(Hi) > 20 and AV < 10,000 Isms" 1 , 
which is the range over which we look for trends. A stricter definition would 
adopt the standard DLA criterion of log A'(Hi) > 20.3. Moreover, at large 
AV our sample will include an increasing fraction of intervening absorbers. 



7 SUMMARY 

We have used a sample of 85 proximate absorbers in the SDSS 
with AV < 10,000 km s _1 and log A'(Hi) > 20 to investigate trends 
of metal line strength with velocity separation. Composite spectra 
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are constructed for 3 AV ranges: AV< 3000 km s~\ 3000 < AV 
< 6000 km s -1 and AV > 6000 km s -1 . The samples are further 
divided according to H i column density and QSO luminosity. Our 
main findings are as follows: 

(i) Metal line EWs are largest in the AV < 3000 km s _I compos- 
ite and smallest at AV > 6000 km s -1 . At intermediate velocities 
(3000 < AV < 6000 km s -1 ) the EWs are between the low and high 
AV composites. We interpret this result as caused by higher metal- 
licities at lower relative velocities. 

(ii) Although the intermediate velocity composite indicates that 
enhanced metallicities might exist out to several thousands of 
km s -1 from the systemic QSO redshift, measurements of Sill A 
1526 in individual spectra indicate that the metallicity is most en- 
hanced at AV < km s _I . 

(iii) The difference between the EWs in the lowest (AV < 3000 
km s _1 ) and highest (AV > 6000 km s -1 ) composites is largest when 
only the high N(Ui) (log Af(Hi) 3> 20.7) are considered. Using the 
Sill A 1526 EW-metallicity relation of Prochaska et al. (2008a), we 
determine metallicities of 1/12 and 1/45 for the low and high AV 
composites respectively, when A'(H i) > 20.7. 

(iv) The absorbers within 3000 km s _I of high luminosity QSOs 
have lower metal line EWs than those within 3000 km s -1 of low 
luminosity QSOs. We speculate that this may be due to truncation 
of star formation in the vicinity of a more luminous AGN. However, 
a low AV and high QSO luminosity does not necessarily lead to 
lower EWs. 

(v) Using the classifications of Herbert-Fort et al. (2006) we 
have also shown that PDLAs are 5 times more likely to be 'metal- 
strong' than intervening DLAs. We present 4 PDLAs whose Sill 
A 1808 is strong enough to be be detected even in the SDSS spec- 
tra which are also classified as MSDLAs, making them excellent 
candidates for follow-up of rare species. 

(vi) No Lyo- emission is detected in the PDLA trough of the 
composite. Although an accurate flux limit is difficult to model 
based on the uncertainty of line emission properties, our simula- 
tions indicate that typical Lyo- luminosities of a few xlO 42 would 
have been detected. 

(vii) From an analysis of the power law slopes of QSOs with 
PDLAs compared to a matched control sample, we place an upper 
limit of E(fi - V) < 0.014 (for an SMC extinction curve) on the 
reddening caused by PDLAs. This is less than inferred from a sim- 
ilar analysis of MS DLAs, and a complementary investigation of 
reddening in proximate Mgll absorbers at z ~ 1.5. 

Our findings are consistent with measurements in a smaller 
sample (16 PDLAs) of high resolution spectra which find that high 
N(H i) PDLAs are more metal-rich than the intervening population 
by a factor of ~ 3 (Ellison et al. 2010). Taken together, these re- 
sults are evidence that galaxies in the vicinity of QSOs not only 
differ from the intervening population but are directly affected by 
the proximity of a QSO. 
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APPENDIX A: PDLAS IN THE SDSS SAMPLE 

In Figure I Al I we present the fits to the Lya absorption of the 85 
PDLAs in our sample. 
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Figure Al. Fits to the Lya absorption in our PDLA sample, the blue vertical dashed line indicating the centre of the Lyo- trough. QSO names are abbreviated 
to Jhhmm+ddmm in order to fit the panels. Full SDSS QSO names are given in Table [T] In cases where 2 PDLAs are present, they are shown in separate 
panels but have been fitted simultaneously due to the significant overlap in velocity space. The grey shaded region shows the uncertainty in the fit. The grey 
dot-dashed line shows the zero level and orange dotted line line the lcr error array. The purple dashed line shows the fit to the continuum. 
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Figure Al. Continued. 
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Figure Al. Continued. 



Proximate DLAs in the SDSS 



J 1257-0 11 


1 z=4.022 N(HI)=2Q.30 

1 /\\ [Y 

\. i / y m 1 

n \ i /j>J lj 


6050 


6100 6150 
Wavelength (A) 




./1419+5923 z=2.247j"N(Hl)=20!95 


M 


Ifjrfk i f$ 
TO i / 






4950 5000 5050 5100 5150 5200 
Wavelength (A) 

J 1 353-0250 z=2.362j N(H1)=2030^ 



3900 3950 _ 4000 
Wavelength (A) 




5250 5300 
Wavelength (A) 



4050 41 00 

Wavelength (A) 



4900 4950 
Wavelength (A) 




4900 4950 
Wavelength (A) 



4600 4650 /700 4750 

Wavelength (A) 



4000 _ 4050 

Wavelength (A) 



J1304+1202 z=^2.928 | N(HI)=20.35 

I 
I 



J 1358+0349 z=2.854,N(HI)=20.40 




4750 4800 
Wavelength (A) 



4650 4700 
Wavelength (A) 



5050 5100 5150 5200 
Wavelength (A) 




4750 4800 
Wavelength (A) 



4900 4950 

Wavelength (A) 



4150 4200 4250 

Wavelength (A) 



30 
20 
10 

oi: 



J1322+5841 z=2.819 | N(HI)=20.50 




4600 4650 

Wavelength (A) 



3950 4000 
Wavelength (A) 



4300 4350 
Wavelength (A) 



Figure Al. Continued. 
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Figure Al. Continued. 



